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ABSTRACT 

This paper contributes to the magnetic bearing 
literature in two distinct areas: high temperature and 
redundant actuation. Design considerations and test 
results are given for the first published combined 538°C 
(1000°F)-high speed rotating test performance of a 
magnetic bearing. Secondly, a significant extension of 
the flux isolation based, redundant actuator control 
algorithm is proposed to eliminate the prior deficiency of 
changing position stiffness after failure. The benefit of 
the novel extension was not experimentally 
demonstrated due to a high active stiffness requirement. 

In addition, test results are given for actuator failure tests 
at 399°C (750°F), 12,500rpm. Finally, simulation results 
are presented confirming the experimental data and 
validating the redundant control algorithm. 

INTRODUCTION 

Large advancements in gas turbine engine efficiency 
mandate continuous operation at temperatures 
substantially higher than the present generation. 
Rotordynamic stability and response present challenges 
in these temperature regimes due to limitations in 
conventionally lubricated bearings and dampers. A host 
of innovative devices ranging from damping seals to dry 
friction dampers, foils bearings and magnetic bearings 
(MB) have been developed as possible solutions. 
Magnetic bearings possess distinct advantage over the 
alternate solutions in their ability to work at high altitudes 
(no air requirements as for foil bearings to adapt to 
various load conditions - take off, landing, maneuvers, 
etc.). The present manuscript provides results from 
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recent milestone testing achieved at NASA GRC for 
combined high speed/high temperature MB operation. In 
addition a novel extension of the flux isolation approach 
to redundant MB actuator control is proposed. 

NOMENCLATURE 

Gy Proportional gain (i, j) = (x, y) 

—p 0S Pos 't' on stiffness matrix 

K Isolated opposing C core position stiffness 

NF 

Kpos Kpos for ur, - failed State 

O , x y control voltage to abc control voltage 

—abc 

distribution matrix 

— a bc Control voltage to a, b, c axes power amplifiers 

LITERATURE REVIEW 

A significant number of research papers have 
recently appeared dedicated to high temperature MB 
development. Such areas as electrical insulation, load 
capacity and thermal growth have formed the focus of 
these efforts. 

Ohsawa, et. al [1] developed and tested a MB at 
399°C (750°F) on an 80kg (1761b), 12,000rpm blower 
rotor. A 1500hr endurance test was successfully 
completed with the high temperature magnetic bearing 
(HTMB). Kondoleon and Kelleher [2] investigated the 
magnetic and structural integrity of magnetic bearing 
laminate materials up to 538°C (1000°F). Satisfactory 
characteristics were obtained for Hiperco 50, Hiperco 
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50HS and Rotelloy 5. Mekhiche, et. al. [3] have 
developed an 8 pole magnetic bearing with capability to 
reach 593°C (1100°F) and operate at 50,000rpm, 
however, only room temperature results were reported. 
Xu, et. al. [4] have developed and tested a 538°C 
(1000°F), 8 pole MB and position sensor system. Only 
static tests were performed at the time of the paper’s 
publication. Scholten [5] performed a simulation study to 
define maneuver load requirements and Hilbner and 
Rosado [6] identified benefits and necessary 
modifications for implementing HTMB’s in gas turbine 
engine applications. 

The two areas to be discussed in the present paper 
are the HTMB magnetic bearing construction and also 
the novel redundant control algorithm used with it for 
magnetic levitation. This algorithm provides fault 
tolerance by maintaining decoupled, linearized magnetic 
suspension, control forces even after failure(s) of coil 
within the actuator or power amplifiers driving current 
through the coils. Previous work on fault tolerant 
magnetic bearing actuators focused on two general 
approaches: flux coupling and flux isolation. The former 
approach as developed by Meeker, Masler, et. al. at the 
University of Virginia [7,8,9] changes fluxes in all failed 
and unfailed poles to provide spatially (x, y) uncoupled 
and linearized control forces, after a failure. The fluxes 
remained coupled between all poles in this approach and 
the maximum flux density limits the load capacity if it 
exceeds the material saturation limit. A derivative 
approach that simultaneously minimizes the Euclidean 
norm of the pole flux densities was developed by Na and 
Palazzolo at Texas A&M [10]. The flux isolation method, 
initially proposed by Storace, Lyons, et. al. [11,12] at 
General Electric, employs a redundant control axis so 
that two of the 3 control axes, in their 12 pole actuators, 
are sufficient to provide the demanded control forces. 
Each axis is composed of two opposing C-Cores that 
have the same, centered actuator, bias flux and equal 
and opposite control fluxes. Both flux isolation and flux 
coupling approaches to fault tolerant MB actuators 
appear in the cited literature. The benefit of the former is 
its simplicity in possessing only 3 current redistribution 
matrices, and its lower hardware requirements (6 power 
amplifiers and no decoupling choke). The drawbacks of 
the former include less robustness for many pole failure 
combinations and lower load capacity. The choice of the 
proper approach for a given application should be based 
on the level of fault tolerance required and consideration 
of the reliability of the power amplifiers and decoupling 
choke, as well as the actuator's coils. 

Significant contributions are provided in the present 
paper for both the HTMB and fail-safe/redundant control 
areas. Specifically the HTMB performance data appears 
to be the first published accounts of an actual rotating 
test (10,000rpm) at 538°C (1000°F), i.e. previous tests 
were performed only statically at temperatures above 


399°C (750°F). It is well known that looseness of parts 
on a spinning shaft may cause rotordynamic instability 
due to internal friction and parametric excitation. The 
success of the rotating test confirmed that the rotor 
lamination fits were effective in preventing anomalous 
vibration at high speed and temperature. It is also 
notable that this was accomplished with a twelve pole 
actuator that could be used with either a flux isolation or 
flux coupling approach to redundant actuation. A novel 
extension of the flux isolation method for redundant 
actuation is also presented, in which the deficiency of 
coupled (x,y) position stiffnesses after failure is 
eliminated. 

REDUNDANT CONTROL ALGORITHM 

A redundant control algorithm was formulated to 
provide stable control during control axes failures. 
Specifically, the purpose of the redundant control 
algorithm is to preserve control forces and effective 
stiffness following a single axis failure. This was 
accomplished by developing current distribution and 
cross-coupled gain matrices for each of the three failure 
cases. 

The 3 axis, isolated flux proposed in [11,12] 
preserves control forces following a single axis failure. 
This method employs the following distribution matrices 
to convert xy control signals into control signals along the 
3 redundant axes (a, b, c), which are in turn routed to 
power amplifiers (Figure 1). 


12 Pole Radial Magnetic Bearing Control Block Diagram 



Figure 1. Feedback Control Block Diagram. 
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These matrices are obtained from a least squares 
minimization of control effort. 


stiffness, K act , and position stiffness, K pos as shown in 
equation (11). 


Labe 



No failure Q . 

—a be 


.577 

.333 

0 

.667 

-.577 

.333 


(1) 

( 2 ) 


K cff = K ac ‘ + KP° S (11) 

Let K cur represent the current stiffness for a single, 
opposing, C core pair, then the active stiffness is given 
by equation (12). 


“a” axis failure 


“b” axis failure 


“c” axis failure 


—abc 


—abc 


—abc 


0 0 
.577 1.0 

- 1.155 0 

.577 


0 

-.577 

1.155 

-.577 

0 


(3) 

(4) 

(5) 


Although control forces are preserved, position 
stiffness related forces are changed by the failure of a 
control axis and are not automatically compensated by 
the above redistribution of currents, i.e. a control law 
change in the PID level must also be activated. The 
following summarizes the change in position stiffness 
with axis failure. Let K represent the position stiffness 
along any control axis (a, b, c) prior to failure. The 
position stiffness matrix for the xy axes has the following 
form; 


Kpos 


No failure 
“a” axis failure 


X pos pos 

^ xx 

is pOS jy- pOS 

_yx ^ yy 

k n p : s =k 


K pos =K 


1.5 0 

0 15 
0.75 -.433 

-.433 1.25 


( 6 ) 

(7) 

( 8 ) 


K act = G p K cur (12) 

In order to maintain the desired effective stiffness, the 
proportional gains after failure G pij are related to their 
counterparts before failure G pi j by: 


where 


^ Y cur _i_ Y^F js 

u pij A ^ ^ pos Jj A posjj 

(13) 

ij~ is 

^cur 

(i,j) =(1,1) = (X, X) 

(14) 

(i. j) = (1, 2) = (x, y) 

(15) 

(i,j) = (2, 1) = (y, x) 

(16) 

(i,j) = (2, 2) = (y, y) 

(17) 


Note that G pij is fully populated whereas G pij - is 

typically diagonal for uncoupled SISO type control. This 
necessitates a proportional cross gain in the PD stage of 
Figure 1 to cancel the cross-coupled position stiffness 
after an axis failure. The effects of changes in the 
position stiffness matrix after failure can also be 
minimized by selecting high controller gains such that the 
active stiffness is the dominant component. Although the 
form of equation (13) is quite simple relative to more 
analytically based algorithms, i.e. H*,, fuzzy logic, sliding 
mode, it correctly compensates for the position stiffness 
changes with failure, and may be easily implemented in 
an industrial PID based controller. Its implementation 
does require alteration to an uncoupled PID controller 
coding due to the xy cross coupled gain and adaptation 
with failure. 


“b” axis failure 
“c” axis failure 


K P os=* 


Kpos=« 


1.5 

0 

.75 

433 


0 

1.5 

. 433 " 

1.25 


(9) 

( 10 ) 


A control law algorithm for maintaining the desired 
effective stiffness before and after failure and canceling 
cross-coupled position stiffness is outlined for PD 
control. Active magnetic bearing properties can be 
represented by effective stiffness and damping. The 
effective stiffness is comprised of the active (control) 


BEARING DESIGN 

The HTMB stator (Figure 2) has 12 poles and is 
composed of 214 laminates, each .35mm (14mil) thick, 
and made of heat-treated Hiperco 50A. The pole width 
is selected to yield the maximum force as determined by 
the pole’s cross sectional area, maximum current 
available and number of turns. The design load capacity 
based on a saturation B of 2.0T and activation of 4 poles 
is 1500N (3371b). This approximation ignores fringing 
and material path reluctance. The rotor lamination’s ID 
and OD are 6.0cm (2.36in) and 7.45cm (2.93in), 
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respectively, and the air gap at the poles is 0.45mm 
(17.7mil). 



Figure 2. Stator Lamination Stack. 

The coils are wound from nickel clad copper and are 
potted with a magnesium oxide compound. Each coil is 
hand wound, cured at 121°C (250°F) for 3 hours and 
retained on its pole by metal pins (Figure 3). The 
resistance of each coil is measured vs. temperature up 
to 815°C (1500°F), as shown by Figure 4. Along with hi- 
potting, this has proven to be a good indicator of coil 
integrity. 




(b) 

Figure 3. HTMB Stator (a) and (b) Rotor. 


EFFECTS OF TEMPERATURE ON COIL RESISTANCE 



0 200 400 600 800 1000 

Temperature (C) 

Figure 4. Resistance vs. Temperature for HTMB 
Coil. 

For a flux isolation (C core) approach to redundant 
magnetic bearings it is essential that the flux in any given 
C core does not leak to other C cores even if 1 or more 
of the C cores are turned off. A 3D FE field analysis was 
performed with only the bias flux density circulating in all 
C-cores. This study confirmed the flux isolation 
characteristic, as demonstrated by Figure 5a (all C-cores 
active) and Figure 5b (5 out of 6 C-cores active). It is 
clear from Figure 5b that there is very little flux leakage 
into the failed C-core from the active C-cores, hence 
demonstrating continued isolation of the C-core fluxes 
after failure. The model in this figure has a smaller L/D 
ratio than the actual HTMB rotor in Figure 3, however the 
isolation results should still be applicable. 
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(b) 


Figure 5. Three Dimensional FEA Prediction of 
Flux Field in (a) Unfailed and (b) Failed HTMB. 


HTMB COMPONENT TESTING (STATIC) 

A test fixture was constructed to measure load 
capacity of the magnetic bearing vs. temperature under 
non-rotating condition (Figure 6). Strain gage load cells 
support the stator and measure the force it exerts on a 
dummy rotor shaft which is centered rigidly with dead 
centers. A cooling jacket is employed to protect the load 
cells by removing heat conducted from coil heaters 


wrapped around the OD of the back iron (Figure 7). The 
assembled MB is encapsulated in high temperature 
insulation to retain heat in the rotor, stator and coil 
components, and current is driven through the coil via 
PWM, servo power amplifiers. Figure 8 shows plots of 
MB force vs. current at 6 selected temperatures in the 
range from ambient to 538°C (1000°F). The force 
reduction with increased temperature results from the 
temperature dependence of the B-H characteristics for 
the stator and rotor lamination stack material. The 
temperatures in this figure are measured between poles, 
on the inside surface of the stator lamination stack, at 
the center of the bearing. It is interesting to note that the 
rotor and coil temperatures were measured to be 530°C 
(986°F) and 690°C (1274°F), respectively, with a stator 
temperature of 543°C (1010°F) and current of 12 amps. 
Figure 9 shows a view of the bearing following the 
component testing. 



Figure 6. Test Fixture for Static Force Temperature 
Test of HTMB. 
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HTMB SYSTEM TEST (DYNAMICS) 

Band type heaters and insulators were wrapped 
around the HTMB stator, which was then rigidly fixed in 
the rotating test rig shown in Figure 10. The shaft weight 
is approximately 114N (25.61b), and its CG is located at 
the HTMB centerline in the axial direction. The shaft 
bearing at the turbine driver end is a duplex pair 
consisting of two ABEC 7, angular contact bearings, 
oriented to provide radial support and an angular 
(moment) release. This bearing set was lubricated with 
Krytox grease for high temperature operation. Eddy 
current position sensors are air-cooled and are 
positioned about 5cm (1.97in) axially offset from the 
HTMB ends. Opposing pairs are employed to 
compensate for thermal expansion of the shaft. The 
auxiliary (backup) bearing is located at the outboard end 
of the shaft, is composed of another ABEC 7, duplex 
pair, is lubricated with Krytox grease and has a 0.175mm 
(6.9mil) radial gap. Currents to the HTMB are supplied 
by six 160V, pulse width modulated, power amplifiers. 
The shaft position signals are routed to a digital control 
unit, which in turn produces control signals for the power 
amplifier. The rig was initially spun to 14,000rpm at room 
temperature with magnetic suspension provided by the 
HTMB. A similar test at 538°C (1,000°F) (stator) and 
with a maximum speed of 10,000rpm was also 
performed. The corresponding coil and turbine end 
bearing temperatures were 677° C (1250°F) and 163°C 
(325°F), respectively. 


Figure 10. Photograph of Rotating Test Rig at 
NASA GRC with Test Bearing. 


FAILSAFE TESTING 

Failsafe tests were conducted to evaluate the 
redundant control algorithm at temperature and speed. 
Each of the three control axes were individually failed 
and the response was recorded using a Nicolet data 
acquisition system. The results indicate that the 


Figure 7. Coil Heater Wrapped Around Stator for 
Temperature Controller. 


Figure 8. Magnetic Bearing Forces vs. Current 
and Temperature for Static Testing. 


Figure 9. Close-up View of HTMB Following 
Static Load High Temperature Test. 
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redundant control algorithm provides stable control 
following any single axis failure event. 

The failsafe tests were conducted at 399°C (750°F) 
and 12,500rpm. A digital switchbox was used to initiate 
a failure along a user-selected control axis. The switch 
action opens the control output circuit to the coil and 
simultaneously sends a TTL to the digital controller, 
which triggers the redundant control algorithm. 

The results shown in figures 11 through 15, indicate 
that stable control is maintained during the single axis 
failure events. Figures 11-13 show the test rig data for 
failures of the a-axis. Figure 11 displays the servo- 
amplifier currents arranged in their respective geometric 
positions. The redistribution of the control currents is 
visible when the switchbox de-energizes both a-axis 
coils. Figure 12 shows the displacements at the rotor 
CG and the outboard catcher bearing, where the failure 
in indicated by the vertical line. Notice the glitch at the 
instant of coil failure. Figure 13 exhibits the rotor orbits, 
also at the CG and outboard catcher bearing. Although 
the control forces are conserved, there is a noticeable 
effect on the orbit geometry. 


PWM CURRENTS 399 X. 12500 ipm 
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ROTOR DISPLACEMENTS: 309 *C. 12500 rpm 
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Figure 12. CG and Outboard Displacements, a- 
Axis Failure 


OUTBOARD ORBIT 



Figure 13. CG and Outboard Orbits, a-Axis 
Failure 


Test results validate the redundant control algorithm 
Figure 11. Coil Currents, a-Axis Failure during all three single axis failure scenarios. It should be 

noted that control is also maintained when the failed 
axes are re-energized. Similar results were obtained 
with both b and c axes failures as seen by figures 14 and 
15. 


7 


Copyright © #### by ASME 


OUTBOARD ORBIT 



stiffness. This is illustrated in the figure by the elements 
that possess two radial divisions, the inner representing 
the stiffness diameter and the outer the mass diameter. 
Both linearized and nonlinear magnetic actuator force 
relations were used for comparison purposes. 

Gravitational and gyroscopic effects were included in the 
model. A classical PD controller with lead-lag 

compensation was designed to be identical to the actual 
digital controller. A lag compensator was used in lieu of 
integral control to eliminate steady-state error. 
Coordinate transformations discussed previously were 
used to output control signals to the actuator's abc 
control axes. Control axis failures were simulated by 
discontinuing current to the failed axis coils and applying 
the appropriate transformation matrix at a pre-selected 
time step. 


Figure 14. CG and Outboard Orbits, b-Axis 
Failure 


OUTBOARD ORBIT 
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Figure 16. FE Model of Test Rotor 


CG ORBIT 
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Figure 15. CG and Outboard Orbits, c-Axis 
Failure 


ANALYTICAL ANALYSES/SIMULATIONS 

An analytical model was developed for theoretical 
analysis and simulation using Matlab in an effort to 
simulate and verify the redundant control theory. A 4 
DOF per node finite element model for the rotor 
assembly combined with a classical controller was 
employed. Several analyses were conducted including 
eigenvalues, steady-state imbalance response, transient 
imbalance and levitation simulations. 

The finite element rotor model, depicted in Figure 16, 
was created using a 1-D mesh of sixteen Euler beam 
elements. Four degrees of freedom per node were 
accomplished by eliminating the axial and torsional 
dynamics. The actuator lamination stack was included in 
the elemental mass but not as additional bending 


Several analyses were conducted including closed- 
loop eigenvalues, steady-state imbalance and transient 
simulations. Eigenvalue analysis closely predicted the 
two experimental critical speeds below 20,000rpm, 
approximately 6,300rpm and 14,200rpm. Although the 
state of imbalance was not experimentally measured 
axis failure simulations were performed with a static 
imbalance eccentricity, u = 2.54e-3mm (0.1 mil). The 
steady-state response to imbalance is given in Figure 17. 
Transient simulations of the response to imbalance 
indicate that stable control is maintained in the event of a 
c-axis failure (Figures 18 and 19). As shown, the failure 
has a negligible effect on the displacements and orbit 
pattern. The cross-coupled stiffness cancellation is not 
activated in this scenario because the active stiffness is 
an order of magnitude higher than the position stiffness. 
In order to illustrate the effects of the cross coupled 
stiffness cancellation, the proportional gains must be 
lowered from their nominal values to reduce the active 
stiffness to the same order of magnitude as the position 
stiffness. 
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Y ca ( mm ) (D Outboard Sensor Amplitude (mm) 


Steady-State Imbalance Response 



ure 17. Steady-State Imbalance Response 


HIGH TEMPERATURE TEST RIG 
Linear, C-Axis Failure, 12,500rpm 
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Figure 18. C-Axis Failure Simulation 


HIGH TEMPERATURE TEST RIG 
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Figure 19. C-Axis Failure Simulation 


SUMMARY/FUTURE WORK 

This paper presents results of development efforts 
for a novel HTMB. The bearing exhibited good load 
capacity, although it is degraded somewhat at higher 
temperatures. Total run time on this bearing above 
371°C (700°F) to this date is 27.5 hours. A novel 
extension of the 3 control axis, isolated flux, redundant 
control algorithm is presented to compensate for position 
stiffness changes that occur after failure. The control 
algorithm was verified through simulations and 
experiments. Dynamic tests at temperature and speed 
indicate stable control is maintained during single event 
axis failures. Although the benefit of the novel extension, 
i.e. compensation for position stiffness change, was not 
experimentally demonstrated due to a high active 
stiffness requirement for these tests, It is felt that it may 
be beneficial for other applications, e.g. a vertical rotor. 
A high load (4450N (10001b) capacity)/high temperature 
538°C (1000°F) magnetic bearing is nearly completed 
and will be tested in the next several months at NASA 
GRC (Figure 20). 



Figure 20. High Load/High Capacity C-Core 
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